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In  various  interpolation  problems  of  experimental  and  pointwise 
computer-constructed  data  the  convexity  and  smoothness  (continuous 
differentiability  to  the  second  order)  of  the  interpolant  is  either 
implied  by  the  nature  of  the  underlying  process  or  expressly  desired. 

It  is  well  known  that  classical  (Lagrange,  Hermite,  etc.)  and  ordinary 
spline  interpolation  procedures  do  not  in  general  produce  convex  inter- 
polants  for  convex  data  apd  more  often  than  not  introduce  spurious 
oscillations  between  data  points.  We  introduce  here  a technique  for 
interpolation  and  display  of  convex  data  without  such  drawbacks.  Its 
basis  is  provided  by  an  explicit  solution  of  the  convex  two-point 
Hermite  interpolation  problem  by  a special  rational  function.  To  meet 
the  smoothness  requirement,  a piecewise  rational  function,  called  an  H- 
Spline,  is  constructed.  The  feasibility  of  the  construction  is  equivalent 
to  the  existence  and  uniqueness  of  a solution  of  a nonlinear  system  of 
equations  which  may  be  brought  into  the  fixed  point  form  x = Fx.  Here  x 
is  an  n-dimensional  vector  whose  dimension  depends  on  the  number  of  data 
points  and  whose  components  represent  linear  functions  of  the  slopes  to  be 
assigned  at  the  internal  data  points.  Section  II  contains  the  basic  defi- 
nitions and  preparatory  remarks.  A proof  for  the  existence  and  uniqueness 
of  an  H-Spline  for  any  finite  number  of  data  points,  based  on  Brouwer's 
Fixed  Point  Theorem  and  the  antitonicity  of  the  operator  F,  is  given  in 
Section  III.  The  general  mathematical  background  for  our  considerations 
is  covered  inl»2»3.  A companion  report**  describes  in  detail  the  actual 
construction  of  H-Splines  on  a computer. 

II.  DEFINITIONS 


We  begin  with  the  following  definitions: 

Definition.  A table  T:  (a^,  y^),  i = 0,  1,  n ♦ 1,  such 

that  a = a0<a1  <...  <an<an+j  = b and  s0<Sj  <. . •<sn_i<sn*  where 
si  = (Xi+1  ” yi)/(ai+i  - a^),  i = 0,  1,  ...,  n,  is  called  a convex  table. | 

Definition.  Let  R*  be  the  class  of  rational  functions  of  order 
2 with  at  most  one  finite  pole.  Given  a strictly  monotone  increasing 
sequence  of  real  numbers  a0,  a j,  ...,  an+i,  an  H-Spline  with  joints 
aj,  i = o,  1,  ...,  n + 1,  is  a function  H(x)  defined  for  a = aQ<x<an+ i 
= b that  satisfies  the  following  three  conditions: 


1 Collatz,  L.,  Functional  Analysis  and  Numerical  Mathematics,  (pp.  350- 
361),  Academic  Press,  New  York  and  London,  1966. 

2 

Ortega,  J.M.  and  Rheinboldt,  W.C.,  Iterative  Solution  of  Nonlinear 
Equations  in  Several  Variables,  (pp.  432-446),  Academic  Press,  1970. 

Rail,  L.  B.,  Nonlinear  Functional  Analysis  and  Applications,  (pp.  16-18), 
Academic  Press,  1971. 

Eger land,  W.  0.  and  Wisniewski,  H.  L.,  "Convex  Interpolation  with  H- 
Splines,"  to  appear. 
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a.  In  each  interval  (a^,  *1+1)  for  i = o,  1,  ...»  n,  H(x) 
is  given  by  some  hjeR*. 

b.  H(x)eC2[a,b],  i.e.,  H(x)  is  twice  continuously 
differentiable  on  [a,b]. 

c.  H* ' (x)>o,  a<j^>,  i.e.,  H(x)  is  strictly  convex  on 

[a,b]. | 

Since  hi(x)  * (Aoi+Aijx  + A2i.x2)/(B0i  + Bux),  ai<x<ai+j, 

A2i#o,  H(x),  as  a plane  curve,  is  represented  by  an  arc  of  an  hyperbola 
(or  a parabola  in  case  Bji  » 0)  between  adjacent  joints.  Hence  the  name 
H-Spline  was  suggested. 


Definition.  Let  i and  r be  positive  numbers.  A function  f(x) 
solves  the  "convex  two-point  Hermite  interpolation  problem"  on  the 
interval  [a,b]if (l)f (a)  = o,  f(b)  * 0,  (2) f * (a)  ■ -l,  f*(b)  » r,  and 
(3)f"(x)>o,  a<x*b.| 

% 

According  to  this  definition,  the  function 


B(x)  = B(x;a,b,Jt,r) 


Ar(x-a) (x-b) 
A(x-a)  + r(b-x) 


solves  the  convex  two-point  Hermite  interpolation  problem  for  arbitrary 
a,  b,  a<b,  and  positive  i and  r.  B(x)  is  the  simplest  H-Spline.  We 
note  that 

B”(X)  .-iiMfciii - (1) 

[A(x-a)  + r(b-x)]'J 


III.  EXISTENCE  AND  UNIQUENESS  PROOF 


With  the  definitions  given  in  Section  II,  the  following  theorem 
holds: 

Theorem.  Given  a convex  table  T:  (a i»yi)»  i * 0,  1,  ..., 

n + 1,  and  endconditions  yl,  y'  ,,  y'<s  , sn<yn+i,  there  exists  a unique 
H-Spline  such  that  H(ai)  » yj,  i « o,  1,  n + 1,  H’(a0)  * y^,  10,1,1 

H'(W‘4i-l 

Proof.  Let  y{,  i * 1,  2.  ...,  n,  be  a sequence  of  slopes  such 
that  yo<s0<yi<S! . . •<*n-l<yn<8n<Jrn+l  and  consider  the  function 
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H(x)  * L^x)  + Bi(x),  x*  ai+J,  i - of  1,  n,  (2) 

where  L.(x)  =j.  + s.(x-a.)  and  B.(x)  = B(x;  a.,a.+1  s.  - y.', 

yi+l  " si)*  “(x)  is  strictly  convex  on  [a,b],  H(a.)’=  y. . H'(a.)  = y', 
% 2 i 

and  HfxleC  ra.bl.  Hfxl  IS  an  H-.^n  1 i no  i f nv\A  An  1 \r  « f a.  - • ^ 


and  H(x)eC  [a,b],  H(x)  is  an  H-Spline  if  and  pnly  if  the  continuity 
requirements 


Bi-lCai)  = BV(a.)>o,  i - 1,  2 

are  satisfied.  Using  (1)  and  setting 


(3) 
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(3)  is  equivalent  to  the  existence  of  a solution  of  the  system 

(1  - xj1/2 

X = = f (x  1 

1 (1  - x2)V2  + Kj  ^ 


X.  = 


(1  - Vi)1/! 


Cl  - Xi  + 1)  */2+  V..,!/! 


= fi ^xi-l, xi-#-i^ ' i=2>**‘»  n-1 


x = 
n 


1 + K x ,V2 
n n-1 


= f (x  ) 
nv  n-1' 


(4) 


}"  !h*  °P®n  cube  V 0<*i<l.  i - 1.  2,  ...n.  This,  in  turn,  is  equiva- 
lent  to  the  existence  of  a fixed  point  in  In  of  the  mapping  F:  I CRn-*Rn 
where  for  x = (Xj,  x2,  ...,  xn)  Fx  is  represented  by  the  columnnvector 


fj(x) 


Fx  ■ I ^j,(x) 


VnW/ 

The  mapping  F has  the  following  properties: 

(P,)  F is  continuous  on  I . 
v 1 n 

<P2>  FInCIn* 

(P-1  F is  antitone  on  I . 
v 3 n 

(P4)  If  v,  wel  and  v<w,  then  v « Fw  and 
w = Fv  imply  v * w. 

(Pi)  and  (P?)  are  obvious,  and  (P3)  follows  from  the  mean  value  theorem 
and  3f. (x)/ax. < o,  i,  j » 1,  ....  n,  xeln<  To  prove  (P4),  we  observe 
that  tne  constants  K*,  i = 1,  n,  can  be  expressed  by  the  components 
of  v and  w as  follows: 


1 v 


K.  . 1 ~ VJ  (i  ~ *i*Al/2  1 ~ "i/1  • vi 

1 vi  l wi-l  / wi  1 vi-l 


1 ” 2,  it*) »n  " 1, 


1 - v 1 - w 

K Jlw"1/2 1 


r1/2 2.  v-1/2 

n-l  w n-1 
n 


The  equality  of  the  product  of  the  left  sides  with  that  of  the  right 
sides  in  (6)  yields,  after  cancellation, 

■ <>» 


a contradiction  unless  * w^,  i = 1,  n,  i.e.,  v = w. 
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We  show  next  that  F has  exactly  one  fixed  point  in  In.  First,  let 
60>o  be  such  that  60<6j  = min  (m2(l+m)"2;  (1  + M )“2;  1/2),  where  m * min 
(*1»  ....  Kr)  and  M * max  (Kx,  ...»  K„).  If  xeln(60):  60<Xi<l  - 60, 
i = 1,  ...,  n,  then  it  is  easy  to  verify  that  the  components  of  F satis- 
fy the  inequalities 


6o<fj,(x)<l  - 6q,  i ■ 1,  ...,  n. 


(8) 


In(6  ) is  compact  and  convex,  F is  continuous  on  In(60),  and»  by  (8), 
FIn(“o)Cln^o)*  Therefore,  we  may  apply  Brouwer's  fixed  point  theorem 
to  conclude  that  F has  at  least  one  fixed  point  x*  in  In(<50).  Further- 
more, by  the  antitonicity  of  F,  the  iterations 


k+1 


Fw' 


k+1  = Fv' 


0,  1... 


(9) 


with  initial  points  (v°)T  = (60,  ...,  6?),  (w0)T  = (l-60,...,l-60) 
define  iterates  such  that  ([1],  [3],  [4]) 


/°<v1<  ...<  vk+1<wk+1<  ...<  w^w0' 


(10) 


The  limits  lim  v = v and  lim  w = w exist,  v^w,  x*  is  contained  in 

](-.!»  )(  + CO 

the  order  interval  <v,  w>,  and,  by  the  continuity  of  F,  we  have  v = Fw 
and  w = Fv  from  (9).  Hence,  in  view  of  (5)  - (P4),  v = w = x*  is  the 
only  fixed  point  of  F in  ^(6^).  Since  the  argument  can  be  repeated 
with  an  arbirtary  positive  6<6q  and  In(6o)Cl  (6),  it  follows  that  F 
has  no  other  fixed  point  in  I besides  x*.  The  unique  H-Spline  corres- 

T n - 

pondence  to  (x*)  = (xj,  x*,  ...,  x*)  is  given  by  H(x)  = H(x)  in 
(2)  with  y!  = x?  s^  + (1  - x?)s^,  i * 1,  ...n.  This  completes  the 
proof  of  the  theorem. | 


NOTE: 

An  outline  for  the  proof  of  the  existence  and  uniqueness  of  ap  H-Spline 
for  a given  convex  table  was  presented  at  the  Twentieth  Conference  of 
Army  Mathematicians  at  the  U.S.  Army  Natick  Laboratories,  Natick, 
Massachusetts,  May  1974. 
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